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Abstract 



A scheme of q- deformation of nonlinear maps is introduced. As a 
specific example, a g-deformation procedure related to the Tsallis q- 
exponential function is applied to the logistic map. Compared to the 
canonical logistic map, the resulting family of q-logistic maps is shown to 
have a wider spectrum of interesting behaviours, including the co-existence 
of attractors - a phenomenon rare in one dimensional maps. 
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1 Introduction 

Emergence of the so-called quantum group structures in certain physical prob- 
lems has led to studies on several g-deformed physical systems [1]. Enthused 
over this, and inspired by the elements of Tsallis statistical mechanics [2], we 
suggest a scheme of g-deformation of nonlinear maps. We then elucidate the 
general features of a g-deformed logistic map related to the Tsallis g-exponential 
function, as a concrete illustration of the scheme of g-deformation of nonlinear 
maps. 

Theory of quantum groups turned the attention of physicists to the rich 
mathematics of g-series, g-special functions, etc., with a history going back to 
the nineteenth century [3]. The g-deformation of any function involves essen- 
tially a modification of it such that in the limit q — > 1 the original function is 
recovered. Thus there exist several g-deformations of the same function intro- 
duced in different mathematical and physical contexts. Here, we are concerned 
mainly with the g-deformation of real numbers and the exponential function. 

Originally, in 1846 Heine deformed a number to a basic number as 
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such that [n] q — ► n when q — ► 1. In 1904 Jackson defined a ^-exponential 
function given by 

oc n 
n=0 

with 

[n],! = [n],[n - 1], • • • [2],[1],, [0],! = 1, (3) 
as the the solution of the ^-differential equation 



dfjx) _ f(x) - f(qx) 
d q x (1 — q)x 



/(*)■ (4) 



It is seen that E q {x) — > exp(a;) in the limit q — ► 1 when the Jackson q- 
differential operator d/d q x also becomes the usual differential operator d/dx. 

The mathematics of quantum groups necessitated a new deformation of num- 
ber as 

q n - q- n 



q-q 1 



(5) 



which also has the required property that in the limit q — > 1, [n] q — > n. The 
associated ^-exponential function is given by the same equation (2) but with 
[n] q defined according to (5). 

In the nonextensive statistical mechanics of Tsallis [2] , a new ^-exponential 
function has been introduced as given by 



e 



9 



(l + (l-q)x) 1 ^\ (6) 



which satisfies the nonlinear equation 

df(x 
dx 



- (mr, (7) 



and has the required limiting behaviour: e q — ► cxp(x) when q — > 1. This 
e q plays a central role in the nonextensive statistical mechanics by replacing 
exp(a;) in certain domains of application; it should be noted that it is natural 
to define a generalized exponential function as in (6) if we consider the relation 

/ T \ N 

e x = lim (1 + -) , (8) 

N ►oo V N J 

and regard 1/N as a continuous parameter. The formalism of nonextensive 
statistical mechanics has found applications in a wide range of physical prob- 
lems [2], including the study of nonlinear maps at the edge of chaos. Here we 
derive another deformation of numbers based on the Tsallis g-exponential func- 
tion defined by (6), and use it to study a g-deformed logistic map as an example 
of the general scheme of g-deformation of nonlinear maps. 
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2 A g-deformation scheme for nonlinear maps 

The series expansion of e x q has been presented in [4] as 

00 r\ n 

< = 1 + £^F< 0) 

n=l 

with 

Q„ = l(g)(2 g -l)(3g-2)---(ng-(n-l)) ) n = 0, 1,2, ■ ■ ■ . (10) 
Let us take 

l-g = e, (11) 

and write 

e% = T e (x) = {l + ex) 1 ' e . (12) 
From (9) the series expansion of T e (x) follows as 

^) = £^ as) 



, n! 

n=0 



with 

T„ = 



1, for n = 0, 

1(1 - e)(l - 2e) • • ■ (1 + (1 - n)e), for n > 1. 



(14) 



Comparison of (2) and (13) suggests that we have here another deformation of 
numbers : 



such that 



Then, we have 



lim [n] £ = n. (16) 



e — >0 



*-§w <17) 

exactly analogous to the expression for E q (x) in (2). 

It is usual to extend the deformation rule for integers, such as in (1) or (5), 
to any X by substituting X for n. For example, the commutation relations of 
the g-deformed su(2) algebra are given by 

JqJ± — J±Jo = ±Jbi J+J- — J-J+ = [2Jo] q , (18) 
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where [2Jo] g is obtained from the definition of [n] q in (5) by replacing n by 2Jo, 
i.e., [2Jo]«j = (q 2Jo — q~ 2Jo )/(q — q^ 1 )- Thus, we take that for any real number 

H- " TTTfn 7 "^)' (19) 

Note that [x] e — ► x when e — ► 0. Further, [0] e = and [l] e = 1, as in the case 
of other deformations in (1) and (5). Here after, we shall denote [x] e simply as 
[x] and take it to be defined by (19) unless stated otherwise. 

The g-deformation scheme for discrete dynamical maps we suggest is the 
following. For the one-dimensional map 



L n+1 



f(x n ) (20) 



the g-dcformcd version is 



x„+i = f([x n ] q ), (21) 

where [x n ] q is, in general, any g-deformed value of i n ; essentially a g-dcformcd 
map is obtained by composing the given map with a basic number deforming 
map. For example, the g-deformed logistic map, corresponding to the definition 
of deformed x as in (19), is : 

x n+1 = a[x n ]{l - [x n ]) = F{x n ) = ^1±^-Z^1 . (22) 

We shall refer to the map in (22) simply as the g-logistic map. When e — ► 
the g-logistic map becomes the usual logistic map. In the following we present 
the interesting properties of the g-logistic map in detail. 



3 The g-logistic map 

Figure 1 displays the q- logistic map function in (22), 

_ a(l + e)x(l-x) 
HX) ~ (l + e(l-aO) 2 ' ( ' 

for a — 4. The usual logistic map corresponds to e = 0. It is clear that as 
e moves away from the map gets more skewed. The map is skewed to the 
right for positive e and to the left for negative e. The lower and upper bounds 
for e are, respectively, —1 and oo where ^(a;) vanishes. Taking the domain of 
the g-logistic map the same as for the logistic map, namely the interval [0, 1], 
it is found that the range of the g-logistic map is the same as for the logistic 
map. This is so because over the interval [0, 1] the map function F(x) has no 
singularity for — 1 < e < oo, and x < [x] < 1 for — 1 < e < and [x] < x for 
e > 0. There is one important qualitative difference between the usual logistic 
map and its deformed version. The deformed map is concave in parts of x- 
space, whereas the usual logistic map is always convex. Note that the form of 
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the g-logistic map is similar to certain maps used to model population dynamics, 
such as the Bellows map : f(x) = rx/(l+x b ), and the Moran-Ricker exponential 
map: f(x) = xcxp{r(l — x)}, where r is the nonlinearity parameter leading 
from periodic behavior to chaos [5] . 




Figure 1: Graph of the g-logistic map function F(x) for a = 4. The solid curve 
is for e = (usual logistic map). The dashed curves skewed increasingly to 
the right in the top panel are for e = 1,5, 25, 100. The dashed curves skewed 
increasingly to the left in the bottom panel are for e = —0.5, —0.8, —0.9, —0.99. 



The fixed points of the map in (22) are given by 

a(l + e)x*(l - x*) 



x = 



{l + e(l-x*)y 



(24) 



•5 



(25) 



This has one solution at x* — 0. For x* ^ the above equation becomes 

_ g(l + e)(l-Q 
(l + e(l-^)) 2 ' 

Substituting 1 — x* = y in this equation gives 

(l + ey) 2 = a(l + e)y, (26) 

or 

eV + (2e-a(l + e))y + l = 0. (27) 

So 



V=^2 {a(l + e)-2e± v / ^(lT7)^4^(TTT)} (28) 
Thus, the fixed points {x*} of the q- logistic map are: 

x* = 0, (29) 

X \ = (i - ^2 + e ) - 2e >) + 2?v^(i + £ ) 2 - Mi + £ )' ( 3 °) 

= (l - ^ {a(l + 6) - 2e}) - ^a 2 (l + - <M1 + 6). (31) 

Figure 2 displays the numerically obtained fixed points of the map for a = 1.5 
for different values of e. These coincide exactly with the analytical expression 
for x* + above. It should be noted that for a given value of a only for some range 
of values of e both x* + and x*_ lie within the range of F(x). 
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Figure 2: Bifurcation diagram of the map with respect to e for a = 1.5 

Figure 3 gives the bifurcation diagram of the map for a = 3 with resepct to 
different values of e. It is clear that the fixed point x* + loses stability for negative 
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Figure 3: Bifurcation diagram of the map with respect to e for a = 3 

e. This can be straight-forwardly understood from the absolute magnitude of 
F'(x*), which crosses the value 1 at e = when a = 3. 

Figure 4 gives the bifurcation diagram of the map for a — 3.5 with respect to 
different values of e. In this case the usual logistic map (e = 0) yields a 4-cycle. 
As e increases there is a reverse bifurcation and the 2-cycle gains stability, 
followed by the fixed point. However, note that the fixed point at x* = 
also gains stability when e is sufficiently high, and we have a co-existence of 
attractors, namely, the fixed point at co-exists, first with the 2-cycle and then 
with the fixed point x* + . 
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Figure 4: Bifurcation diagram of the map with respect to e for a = 3.5 
Similar features, namely reverse bifurcations in e-space, and the co-existence 
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of the fixed point x* = with other dynamical behaviour at high e, are observed 
for larger values of a as well (see Figs. 5-9). 
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Figure 5: Bifurcation diagram of the map with respect to e for a = 3.6 
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Figure 6: Bifurcation diagram of the map with respect to e for a = 3.7 



4 Co-existence of attractors 

It is clearly evident from the above that the g-logistic map offers a rare example 
of multiple attractors in a one-dimensional smooth unimodal system [6]. The 
fixed point at x* — co-exists with other kinds of dynamical behaviour when 
e is sufficiently high. The basin of attraction for x* — 0, as reflected in the 
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Figure 7: Bifurcation diagram of the map with respect to e for a = 3.8 

fraction of initial conditions which are attracted to 0, monotonically increases 
with e. Figures 10-12 show the basins of attraction for the fixed point x* = for 
three different values of a. It is clear that this fixed point is a global attractor 
at sufficiently high e. In Figure 12 for a — 4 clearly there is a sharp transition to 
a global attractor (where all initial conditions lead to x* = 0) at around e ~ 3. 
So, the fully chaotic logistic map under deformation with positive e can yield 
stable fixed points. 

The co-existence of attractors is also evident from the bifurcation diagrams 
for the map with respect to a for different values of e displayed in Figs. 13-15. 

Figures 16-18 display the Lyapunov exponents, obtained from trajectories 
arising from different initial conditions. Fig. 16 displays the exponents with 
respect to e for a = 3.6, and it shows two branches of Lyapunov exponents 
arising from different initial conditions, after e ~ 2.5. Both these co-existing 
exponents are below 0, indicating the co-existence of two regular dynamical 
behaviours. This is bourne out in Fig. 5 - from where it is clearly evident that 
after e <~ 2.5 the fixed point x* = co-exists with a 2-cycle. 

In Fig. 17 one observes, for certain e, both positive and negative Lya- 
punov exponents are obtained for different initial conditions. This signals a 
co-existence of chaos and regular dynamics, as is clearly bourne out by the 
bifurcation diagram in Fig. 6. 

From the Lyapunov exponents shown in Fig. 18 and the corresponding bi- 
furcation sequence in Fig. 9, it is clear that for a = 4, we have a global chaotic 
attractor extending over the entire interval for the usual logistic map (e = 0), 
but obtain a global attractor at x* — for e > 3. So g-deformation of a chaotic 
map can lead to stable fixed points. 

The co-existence of attractors for the g-logistic map can be understood as 
follows. Consider, for example, the case of a = 3.7and e = 7 where a 2-cycle 
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Figure 8: Bifurcation diagram of the map with respect to e for a = 3.9 



co-exists with the fixed point x* = 0. Figure 19 shows the form of the dynamical 
evolution function F{x) with respect to x for these values of a and e. The figure 
also displays the F(x) — x line and the F(x) = x*_ line. The intersection of 
the former with the F(x) curve yields the fixed points x* = 0, x*_ and x* + . The 
intersection of the latter with the F(x) curve yields the two pre-images for the 
fixed point x*_. One pre-image is simply x*_ (i.e., it is also on the 45° line) as it 
is a fixed point solution. The other pre-image is given by : 



_ x (2ee'x* + a') + y / (2ee'x*_ + a') 2 - 4(e 2 x* + a')e ,2 x*_ 

F {X - } ~ 2(e 2 x*_ +a') l42: 



+ e and a' = ae' . This marks the beginning of the interval 
[F^ 1 (x*_ ), 1] which maps on to the interval [0, x* ] in the consequent iteration. 
These two intervals are dynamically connected and all points in it are attracted 
to the fixed point x* = 0. Notice however that the interval [x* , F^ 1 (x*_)\ is 
never mapped to the other two. In fact it always maps on to itself. Basically the 
F(x) in the range [x* , F^ 1 {x*_ )] will lie in the interval [x*_,F ma ,x], where F max 
is the maximum of the map F(x). Whenever F max < F~ 1 (x*_), the middle 
segment will map onto itself. When this happens the dynamics in this segment 
will be distinct from that in the other two, and we will have a co-existence of 
dynamical attractors. Figure 19 also shows the numerically obtained basin of 
attraction for the fixed point x* = as a black bar on the x-axis, and clearly it 
falls exactly in the intervals outlined above. 



10 




-10 1 2 3 4 5 

epsilon 



Figure 9: Bifurcation diagram of the map with respect to e for a = 4 

5 Conclusions 

To summarize, we have introduced a scheme of g-deformation of nonlinear maps, 
enthused over studies on several q-deformed physical systems related to quan- 
tum group structures and inspired by the mathematical basis of Tsallis statis- 
tical mechanics. We characterized the scheme in detail with reference to the 
logistic map. The resulting family of q-logistic maps has been found to have an 
interestingly wide spectrum of behaviours, compared to the usual logistic map. 
In particular one observed co-existence of attractors - a phenomenon rare in 
one dimensional maps. 

One can also study the g-deformations of maps following other deformation 
schemes, like based on (1) or (5). We feel that the study of such families of de- 
formed maps should be profitable for analytical modeling of several phenomena, 
as one parameter, namely the deformation parameter, can be used to fit a large 
range of functional forms, as evident from Fig. 1. For instance, it is seen that the 
experimentally constructed one dimensional map for the Belousov-Zhabotinskii 
reaction in a stirred chemical reactor [7] has a striking similarity to the q-logistic 
map in Fig. 1 for negative epsilon value. Clearly there are many interesting ar- 
eas of research for further exploration. One can study the g-deformations of 
the various other nonlinear maps (one or higher dimensional) and coupled map 
systems. Further, in higher dimensional cases and coupled map systems one can 
experiment with different deformation parameters (e) for different variables and 
for different maps coupled in the system. In conclusion, g-deformation allows us 
to generate interesting families of maps and has potential utility in constructing 
iterated functional schemes to model low dimensional dynamical phenomena. 
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Figure 10: Fraction of initial conditions attracted to the fixed point x* = 0, for 
varying values of e for a = 3.5. 
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Figure 11: Fraction of initial conditions attracted to the fixed point x* = 0, for 
varying values of e for a = 3.7. 
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Figure 12: Fraction of initial conditions attracted to the fixed point x* = 0, for 
varying values of e for a — A. 




Figure 13: Bifurcation diagram of the map with respect to a for e = 1 
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Figure 14: Bifurcation diagram of the map with respect to a for e = 3 
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Figure 15: Bifurcation diagram of the map with respect to a for e = 5 
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Figure 16: Lyapunov exponents with respect to e for a = 3.6. 
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Figure 17: Lyapunov exponents with respect to e for a = 3.7. 
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Figure 18: Lyapunov exponents with respect to e for a = 4. 




Figure 19: The form of the dynamical evolution function F(x) with respect to 
x for a = 3.7 and e = 7. The dashed lines are: the F(x) = x line (45° line) 
and the F(x) = x*_ line. The figure also shows the numerically obtained basin 
of attraction for the fixed point x* = as a black bar on the x-axis. 
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